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1.  Introduction 


Given  the  type  of  a  threat  and  its  attack  direction  the  vulnerability  of  a  target  system 
(e.g.  of  a  tank)  is  computed  in  two  steps.  First,  one  computes  estimates  of  the  effect  of  the 
threat  for  a  number  of  attack  points.  Second,  these  effects  are  integrated  to  obtain  the 
overall  effect.  The  evaluation  of  the  integrals  is  done  by  a  summation  of  the  local  effects 
and  the  result  of  the  integration  measures  the  overall  "vulnerability"  of  the  target. 
Depending  on  the  definitions  of  the  local  effects  and  the  form  of  the  integrals  the  computed 
result  might  be  a  "vulnerable  area",  a  "probability  of  kill"  or  another  similar  averaged 
vulnerability  indicator.  The  accuracy  of  the  result  depends  on  the  following  three  parts  of 
the  computing  process:  the  accuracy  of  the  local  vulnerability  model,  the  detail  and  accuracy 
of  the  target  description  and  the  accuracy  of  the  numerical  evaluation  of  the  integrals.  In  a 
reasonable  procedure,  these  three  sources  of  inaccuracies  should  be  properly  balanced.  'I he 
present  report  provides  criterions  and  error  estimates  that  help  one  to  achieve  such  a 
balance.  For  instance,  using  the  results  of  this  report  one  can  determine  if  improvement  of 
any  part  of  the  calculations  is  indicated,  or  whether  a  contemplated  improvement  will  have 
a  noticeable  effect  on  the  overall  accuracy  of  the  computation. 

Usually  the  numerical  evaluation  of  the  integrals  is  done  by  defining  a  computing  grid 
with  square  cells  in  a  reference  plane  normal  to  the  direction  of  the  threat,  computing  the 
integrands  at  a  random  point  in  each  cell  and  approximating  the  surface  integrals  by 
appropriately  weighted  sums  over  the  computing  cells.  The  discretization  error  will 
decrease  for  smaller  cell  sizes  and,  therefore  the  cell  sizes  should  be  small.  On  the  other 
hand,  the  detail  of  the  target  is  finite,  and  sometimes  the  target  description  is  quite  coarse. 
This  means,  that  for  a  given  target  there  is  a  bound  for  the  cell  sizes  in  the  sense  that  any 
calculations  with  smaller  cells  do  not  change  the  value  of  the  integral,  because  in  a  very  fine 
net  the  integrand  is  constant  within  most  of  the  cells.  Consequently,  there  is  also  a  bound 
for  the  accuracy  of  the  vulnerability  indicators  which  cannot  be  improved  by  refined 
numerical  integration.  That  bound  depends  on  the  detail  of  the  target  description,  on  the 
spacial  variation  of  the  response  of  the  various  components  of  the  target  to  the  threat,  and 
on  the  accuracy  of  the  local  vulnerability  models.  Our  task  is  to  find  this  bound  in  terms  of 
easily  obtainable  input  information.  The  largest  cell  size  that  produces  the  integral  value 
with  the  achievable  accuracy  can  be  considered  as  an  optimal  cell  size.  Also  useful  are 
accuracy  estimates  for  computations  with  cell  sizes  that  are  larger  than  the  optimum.  Then 
the  accuracy  can  be  improved  by  sampling  at  more  than  one  point  from  each  cell,  and  one 
would  like  to  have  an  estimate  of  the  proper  number  of  samples. 

We  illustrate  some  aspects  of  the  problem  using  a  one-dimensional  example  where  the 
problem  is  easily  tractable.  Let  y=h(x)  be  a  piecewise  constant  function,  and  let  fj  be  the 
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sum  of  lengths  of  segments  where  y-hs  =  constant .  Then  (see  Figure  1) 

b  c 

I  =  fh(x)  dx  =  2/j  hi  ■  (1-1) 

a  j-  1 

Let  the  values  hj  have  uncertainties  which  are  expressed  by  spreads  6 h}.  The  corresponding 
spread  of  the  integral  is 

87  =  £/,8*,  •  (!  .2) 

We  compute  the  integral  numerically  by  subdividing  the  segment  b-a  into  n  equal 
elementary  segments  and  sampling  the  integrand  at  k  random  points  in  each  elementary 
segment.  The  average  distance  between  the  sampling  points  is 


and  the  integral  is  approximated  by 

'  -  ■  a-") 

where  xtl  are  the  coordinates  of  the  sampling  points  in  the  * -th  segment.  We  bound  the 
difference  between  7  and  1  by  considering  points  of  discontinuity  of  h(;r).  Let  {fc}  be  the  set 
of  all  those  elementary  segments  which  contain  at  least  one  discontinuity,  and  let  A  hj  be  ‘he 
range  of  h(x)  within  the  elementary  segment  j.  Then 

|  /  =  ^AA,  .  (1.5) 

"  Mb) 

The  sample  size  parameter  k  does  not  enter  the  expression  (1.5)  for  the  bound  of  the 
difference,  because  the  xlt  are  chosen  at  random  in  each  elementary  segment. 

Next  we  derive  three  estimates  for  the  standard  deviation  o(7)  of  the  approximation  7 
from  the  true  value  I.  Each  estimate  is  based  on  different  assumptions  about  the  integrand 
and  the  arrangement  of  the  sampling  points.  We  show  in  the  Appendix  that  for  a  random 
arrangement  of  nk  sampling  points  the  standard  deviation  is  bounded  by 

«(0  s  yjf-  [  (K..-W-KJ  ]“’  *  j  .  (1.6) 

where  h  is  the  average  value  of  h(x).  This  estimate  only  requires  that  the  integral  /  exists 
and  no  other  properties  of  h(x)  are  used. 

A  better  estimate  is  obtained  if  one  takes  into  account  the  piecewise-constant  property 
of  the  integrand.  Then  the  sum  (1.4)  can  be  rearranged  as  follows 
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(1.7) 


I  =  ( b~a )  hj  , 

j- 1 

where  ay  is  the  fraction  of  sampling  points  that  fall  within  the  segment  fr  The  expected 
value  of  dj  is  fj(b-a )  and  the  standard  deviation  of  a}  is  (Ref.  1,  p.  168) 

CT(ay)  =  \aj  (1-fly)  J1'2  •  (!-*) 


Hence 


or(7)  =  (b-a) 

^oia,)2  ( hj  -  h)2 

1/2 

_  1 

VnT 

i/y(l-/y)(Ay-A)2 

. 

j’1 

1/2 


(1.9) 


For  the  third  formula  we  assume  that  the  integrand  is  piecewise  constant  and  that  the 
n  elementary  segments  are  equal.  In  each  elementary  segment  the  k  sampling  points  are 
random  and  we  may  use  eq.  (1.6)  to  estimate  the  standard  deviation  of  integration  over  the 
segment.  Let  lt  be  the  integral  of  h( x)  over  the  j-th  elementary  segment,  h,  be  the  average 
value  of  h(x)  within  the  segment,  and  7j  =  (7>-a)/i,/n.  Then,  in  analogy  to  eq.  (1.6) 


<>(/,)  -■  i 

Hence  an  estimate  of  the  standard  deviation  of  7  is 

a(T)  = 


1/2 


b-a  1  AI 

TVk  2  4*‘ 


2  «■</■.)  r = ±=t  M  i  w  ] 

[m  nV*  2 


1/2 


(1.10) 


(1.11) 


If  n  is  sufficiently  large  then  the  estimate  (1.11)  is  better  (smaller)  than  the  estimate  (1.9). 
This  is  so  because  in  the  derivation  of  cq.  (1.1!)  we  have  used  more  information  about  the 
computing  process. 

The  number  of  sampling  points  should  be  sufficiently  large  so  that  the  error  of  the 
numerical  integration  is  smaller  than  the  intrinsic  error  of  the  integral.  Therefore,  one 
should  chose  n  such  that 


D  <  6/  ,  (1.12) 

or  chose  n  and  k  such  that,  using  one  of  the  estimates  for  <7, 

a(f)  <  67  ,  (1.13) 

or  less  than  a  fraction  of  87.  The  second  condition  (1.13)  is  less  demanding,  because  it  uses 
an  estimate  of  the  standard  deviation  of  the  integral,  whereas  the  condition  (1.12)  is  based 
on  a  bound  of  the  integration  error. 

A  function  step  with  y^hj  and  the  width  tj  contributes  ( hj-h0)tj  to  the  value  of  the 
integral,  whereby  h0  is  the  average  value  of  y  within  the  segment  tj  if  the  step  is  left  out.  It 
is  not  reasonable  to  include  the  step  into  the  definition  of  the  function  if 
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|  (hj-hq)  t,  I  «  hi 


(114) 


or  if 

|  (hj-h0)  tj  |  «  o(T)  .  (1.15) 

In  the  example  shown  in  Figure  1  the  step  of  the  integrand  with  the  width  /3  contributes  to 
the  integral  less  than  hi  (half  of  the  shaded  area).  Therefore,  this  step  can  be  deleted 
without  a  noticeable  effect  on  the  integral  value.  This  criterion  for  a  deletion  or  inclusion 
of  a  step  into  the  integrand  depends  only  cn  the  accuracy  of  the  integrand  and  is 
independent  of  the  integration  method.  An  example  of  a  border  case  is  the  step  with  the 
width  /2.  Its  contribution  to  the  integral  is  approximately  equal  to  the  spread  hi  and, 
therefore,  it  should  be  kept  in  the  integrand. 

To  capture  all  relevant  parts  of  the  integrand,  the  distance  between  sampling  points 
should  be  less  than  the  width  /  of  the  smallest  relevant  step  of  h(x).  In  terms  of  the 
average  distance  (1.3)  the  condition  is 

(M6) 

n>c  l 

More  rigorous  is  the  condition  in  terms  of  the  maximum  distance  between  sampling  points: 

(M7) 

n  l 

The  conditions  (1.16)  or  (1.17)  for  the  number  of  sampling  points  are  based  on  the  detail  of 
the  function  h( x),  whereas  the  conditions  (1.12)  and  (1.13)  are  based  on  the  accuracy  of  the 
function.  Hence,  each  set  of  conditions  address  different  aspects  of  the  problem  and 
supplement  each  other. 

In  the  remaining  part  of  this  report  we  apply  the  outlined  considerations  to  the  surface 
integrals  of  vulnerability  computations.  To  simplify  the  analysis  we  assume  that  the  number 
of  computing  cells  or  sampling  points  is  large  (at  least  of  the  order  of  100)  since  this  is 
usual  in  applications.  The  integral  formulation  of  vulnerability  indicators  is  provided  in 
Section  2,  where  we  also  discuss  the  inaccuracies  of  the  vulnerability  integrals  in  terms  of 
the  inaccuracies  of  the  predicted  local  effects.  From  that  discussion,  we  establish  criterions 
for  the  inclusion  of  (small)  additional  components  into  a  given  target  description.  Formulas 
for  the  numerical  evaluation  of  the  integrals  by  sampling  are  given  in  Section  3.  The 
section  also  contains  accuracy  estimates  for  the  integrals  based  on  the  number  of  sampling 
points.  In  Section  4  we  discuss  the  accuracy  of  the  numerical  evaluation  of  the  integrals  in 
terms  of  the  properties  of  the  integrand.  Conditions  are  derived  for  the  size  of  the 
integration  cells  and  for  the  number  of  samples  which  insure  that  the  target  geometry  is 
properly  taken  into  account.  The  accuracy  estimates  of  Section  3  are  specialized  in 
Section  5,  where  we  take  into  account  the  piecewise-constant  property  of  the  integrand.  A 
further  refinement  is  provided  in  Section  6,  where  also  the  cellular  arrangement  of  the 
sampling  points  is  exploited.  In  Section  7  we  illustrate  the  application  of  the  results  by 
presenting  an  example  where  the  proper  cell  sizes  and/or  number  of  sampling  points  are 
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estimated.  A  summary  of  the  results  and  conclusions  are  given  in  Section  8. 

2.  Integral  Formulation  of  Vulnerability  Indicators 

We  consider  the  vulnerability  of  a  target  system,  such  as  a  tank  with  respect  to  a 
threat  from  a  particular  direction.  To  the  threat,  the  target  presents  an  area  (a  silhouette) 
within  which  the  threat  might  have  an  effect.  The  silhouette  is  the  projection  of  the  outline 
of  the  target  onto  a  plane  normal  to  the  direction  of  the  threat.  We  call  this  plane  the 
reference  plane.  In  the  vicinity  of  the  target  the  threat  is  assumed  to  travel  along  a  straight 
path  called  a  shotline.  Each  shotline  is  norma]  to  the  reference  plane  and  completely 
described  by  the  coordinates  x  and  y  of  its  intersection  with  the  plane.  If  a  shotline  of  an 
attacking  threat  falls  within  the  silhouette  of  the  target  then  the  threat  might  have  a  non¬ 
zero  effect  which  we  assume  can  be  calculated  and  which  we  denote  by  v(x,v).  We  further 
assume  that  the  effect  function  v(x,y)  is  dimensionless  and  has  a  value  between  zero  and 
one.  The  overall  vulnerability  of  the  target  to  a  threat  from  a  direction  normal  to  the 
reference  plane  is  obtained  by  integrating  the  effect  v(x,y)  over  the  reference  plane.  The 
effect  may  be,  for  instance  the  probability  of  kill  or  the  reduction  of  combat  effectiveness, 
and  the  integral  over  v(x,.v)  can  be  either  simple  or  weighted  (the  latter  if  an  average  of  v  is 
computed).  The  threat  might  also  have  a  distribution  in  the  reference  plane,  specified  h  r 
instance  by  an  aim  point  and  a  dispersion.  Then  the  vulnerability  integral  is  a  weighted 
average  of  v(x,y)  with  the  distribution  density  of  the  threat  as  weight. 

We  define  two  generic  vulnerability  integrals  which  represent  two  types  of 
vulnerability  indicators.  Let  E  be  the  reference  plane,  S  be  the  area  of  the  silhouette  in  E 
and  let  v(x,y)  be  the  probability  of  kill.  We  call  the  first  integral  the  vulnerable  area  and 
define  it  by 


Vs  =  f  v(x,y)  dxdy  .  (2.1) 

5 

The  dimension  of  V s  is  m2. 

The  second  vulnerability  integral  represents  an  average  of  the  effect  function  v(x,y) 
over  an  averaging  area.  That  area  usually  is  either  the  silhouette  of  the  target  or  a 
dispersion  area  of  the  threat.  We  assume  that  the  averaging  area  and  the  averaging  weight 
are  specified  by  a  weight  density  function  g(x, y)  that  has  the  dimension  1/m2  and  is 
normalized  by 


f  g(x,y)  dx  dy  =  1  (2.2) 

E 

We  define  the  average  vulnerability  V A  by 

V a  ~  J"#(-v'.v)  v(x,y)  dx  dy  .  (2.3) 

.4 
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where  A  is  the  area  defined  by  g(x,y)>0.  We  introduce  a  dimensionless  weight  function 
■y(x,y)  by  the  definition 

"i(x,y)  =g( x,y)  •  A  .  (2.4) 

Then  the  average  vulnerability  is 

vA  =  t  / y(x>y)  v(x*y) dx  dy  •  (2-5) 

A  A 

We  assume  that  the  function  g(x,y)  is  exactly  known,  differentiable  in  A  and  has  a 
bounded  gradient  in  A.  About  v(x,_y)  we  assume  that  it  might  contain  unknown  errors  with 
absolute  values  less  than  a  given  error  bound  function  6v(x,>0.  For  simplicity,  we  assume  a 
symmetric  error  bound  with  the  property 

Sv(x,y)  <  min  {  v(x,y),  l-v(x,y) }  ,  (2.6) 

so  that 

0  s  v(x,y)-6v(x,;y)  s  v(x,y)  +  8v(x,y)  <  1  .  (2.7) 

The  intrinsic  error  of  Vs  due  to  the  inaccuracies  of  v(x,y)  is  bounded  by 

8VS  -  f  5v(x,y)  dx  dy  .  (2.8) 

5 

Next,  we  subdivide  the  silhouette  S  in  sub-areas  Fj  such  that  within  each  Fj  the  effect 
function  is  approximately  constant.  Often  these  areas  will  coincide  with  presented  areas  of 
target  components,  particularly  in  so-called  compartment  model  calculations.  In  general 
they  will  be  extensions,  combinations  or  parts  of  component  projections.  Let  the  Ft  be 
called  nodal  areas  and  the  corresponding  constant  values  vFj  of  the  effect  function  be  called 
nodal  values  of  v(x,y).  Formally,  we  define  the  F}  as  follows.  Fjt  j=l,...,c  is  a  set  of 

areas  such  that  if  j*i,  (J  FjQS  and  |  v(x,>)-vf;.|<Avf;/2  for  (x,y)  in  Fp 

t 

whereby  vFj > 0  is  the  average  value  of  v(x,y)  in  Ft,  and  the  are  arbitrary  but  fixed 
constants.  We  supplement  this  set  with  a  nodal  area  F0  which  encloses  all  those  parts  of  the 
silhouette  and  reference  plane  where  vFO=0.  In  terms  of  the  nodal  areas  the  vulnerable 
area  is 

v,  =  2/  dx  dy =  2  Fi  ^ Fj  ■  (2-9) 

jmiFj  l 

The  error  bound  (2.8)  can  be  further  bounded  using  the  FJ  by 

6V,  ^  Ds  =  j^Fj  SvFjm„  ,  (2.10) 

y-i 

where  8vfymil  is  the  maximum  value  of  Sv(x,y)  in  the  nodal  area  F r 
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The  nodal  value  representation  (2.9)  of  the  vulnerability  integral  is  not  merely  an 
artifact.  Such  a  representation  is  natural  for  compartment  model  calculations  where  v{x,y) 
is  approximately  constant  within  the  projections  of  component  silhouettes  onto  E.  Also, 
recent  high  resolution  calculations  that  include  spall  modeling  indicate  that  v(x,y)  likely  in 
general  assumes  only  a  finite  number  of  values.  Examples  and  an  explanation  of  this 
property  of  v(x,y)  are  given  in  Reference  2,  page  35. 


Eqs.  (2.8)  or  (2.10)  can  be  used  to  determine  whether  a  particular  component  should 
be  added  to  an  existing  target  description.  Let  the  effective  presented  area  of  a  component 
be  that  part  of  £  in  which  the  threat  has  an  effect  on  the  component.  In  a  compartment 
model,  the  effective  presented  area  is  the  projection  of  the  silhouette  of  the  component  onto 
£.  In  a  general  model,  particularly  if  spall  effects  are  considered,  the  effective  presented 
area  will  be  larger  than  the  silhouette.  Let  Tj  be  the  effective  presented  area  of  a 
component,  vTJ  be  the  average  value  of  v(x,y)  within  T]t  v0j  be  the  average  within  T} 
before  the  component  is  added,  and  hvTJmtx  be  the  maximum  of  8v  within  Tr  Then  a 
reasonable  criterion  for  the  inclusion  of  the  component  into  the  target  description  is 


or 


T l  ‘  (y tj  v0J)  I  >  2 
T]  '  8vr;mtlI  >  —  Ds  . 


(2.11) 


The  second  condition  in  (2.11)  is  included  to  cover  cases  where  the  addition  of  the  new 
component  does  not  significantly  change  the  value  of  Vs,  but  nevertheless  increases  the 
.curacy  of  the  result.  The  exclusion  of  such  a  component  would  falsely  suggest  a  too 
accurat;  result. 


Corresponding  formulas  for  the  average  vulnerability  VA  are  as  follows. 

8V,  =  Y  fy(x,y)  Bv(x,y)  dx  dy  ,  (2.12) 

A  A 

=  .  (2.i3) 

A  7=0 

8V,  *  Da  =  i  2  (£,n  -4)  ■  bt  8v  ]FJm,x  .  (2.14) 

A  7=0 

The  summation  in  eqs.  (2.13)  and  (2.14)  starts  with  j— 0  because  the  nodal  area  with  vf;  =  0 
must  be  included  in  the  sum.  The  total  number  of  nodal  areas  intersecting  with  A,  e  +  l,  in 
general  will  be  less  than  the  number  c  of  nodal  areas  intersecting  with  the  silhouette.  The 
criterion  for  the  inclusion  of  a  component  with  the  effective  presented  area  7V  in  the  target 
description  is 
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Tj  •  J  •  {  MT}-(yv)0J 
and  Tj  ■  j  ■  fr  8v  ]TJmtx 

Parts  of  nodal  areas  which  satisfy  the  conditions  (2.11)  or  (2.15),  respectively,  we  call 
relevant.  One  can  use  these  criterions  to  simplify  target  descriptions  so  that  all  parts  of 
nodal  areas  are  relevant.  In  particular,  the  criterions  can  help  to  simplify  complicated 
boundaries  of  nodal  areas.  We  note  that  in  certain  problems  there  can  be  other  criterions 
which  justify  the  inclusion  of  small  nodal  areas  in  a  target  description.  For  instance,  if  the 
target  contains  a  large  number  of  irrelevant  components  which  cannot  be  easily  combined 
into  a  few  components  with  simple  geometry,  and  whose  combined  contribution  is  relevant, 
then  one  might  treat  all  components  as  relevant.  However,  as  we  shall  see  later,  if  the 
relevent  details  are  very  small,  then  one  needs  a  very  long  computing  time  to  capture  their 
small  contributions.  Therefore,  it  is  advisable  to  keep  the  target  description  simple, 
smoothing  out  complicated  component  boundaries  and  combining  small  components  into 
larger  and  simpler  ones.  Because  only  approximate  estimates  of  Ds,  DA  and  of  the  other 
terms  in  the  criterions  (2.11)  and  (2.15)  are  needed,  the  estimates  can  be  obtained  by  visual 
inspection  of  temporary  results  with  few  and  simple  components. 

3.  Monte  Carlo  Computation  of  Vulnerability  Integrals 

The  computation  of  a  value  of  the  effect  function  v(x,y)  typically  involves  tracing  the 
shotline  with  the  coordinates  (jc,y)  through  the  target  and  evaluating  the  damage  models  for 
those  target  components  which  are  encountered  by  the  shotline  or  otherwise  affected,  e.g. 
by  spall.  The  results  from  the  damage  models  are  aggregated  by  special  algorithms  to 
obtain  the  effect  v(x,y)  of  the  threat  on  the  target  system.  Because  v(jr,y)  can  only  be 
obtained  pointwise,  the  vulnerability  integrals  are  computed  by  discrete  approximations. 
Usually  this  is  done  by  defining  in  the  reference  plane  a  square  grid,  obtaining  a  value  of 
v(x,y)  for  a  randomly  chosen  point  within  each  cell  of  the  grid,  and  adding  the  values  with 
appropriate  weights  according  to  eq.  (2.1)  or  (2.5),  respectively. 

For  a  general  grid  consisting  of  identical  cells,  the  algorithm  can  be  described  as 
follows.  Let  C  [mJ]  be  the  area  of  a  cell  in  the  computing  grid  and  n  be  the  number  of  cells 
within  the  silhouette.  Let  (x(,y,)  be  a  point  in  the  i-th  computing  cell.  Then  the 
approximation  of  the  vulnerable  area  integral  is 

V*  =  C  2V(W()  =  -  £ ''(*,->,)  -  (3.1) 

i-i  "  /-t 

where  the  summation  is  over  the  n  cells  that  intersect  with  the  target  silhouette.  The 
approximation  of  the  average  vulnerability  integral  is 


>iD‘ 


>iD*  • 


(2.15) 
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(3.2) 


VA  =  ~r  2^-w,)  vCjri»>'i)  =  —  2"^)  v(*,.:y,)  . 

A  i-1  m  i-1 

where  m  is  the  number  of  cells  that  intersect  with  A. 


A  simple  generalization  of  the  algorithms  is  to  take  more  than  one  sample  of  the 
integrand  from  each  computing  cell.  Let  k  be  the  number  of  samples  from  each  cell,  let 
(x,t,ya)  be  points  within  the  i-th  cell  and  let  vu  =  v(xil,yll).  Then  the  approximate  vulnerable 
area  integral  is 


1  n  k  n  k 

2*- 


i-l  /-I 


/- 1  /-i 


(3.3) 


and  the  approximate  average  vulnerability  integral  is 

.  m  k  .  m  k 


•i  /-i 


mk 


V.7 


i-1 /-l 


(3.4) 


Equations  (3.3)  and  (3.4)  include  as  special  cases  eqs.  (3.1)  and  (3.2). 


Estimates  of  the  standard  deviation  of  a  numerically  computed  integral  from  its  true 
value  are  derived  in  the  Appendix.  The  estimated  standard  deviation  of  the  approximation 
Vs  to  the  vulnerable  area  integral  Vs  is 


<r(v»)  s  v* 


s  j_  _  j  Vs 
Vrt  2  2  \Jk  ’ 


(3.5) 


where  s~\/C-(S/n)V2  is  the  side  length  of  a  computing  cell  if  the  grid  consists  of  square 
cells.  An  estimate  of  the  standard  deviation  of  the  approximation  VA  to  the  average 
vulnerability  VA  is 


°(Ya) 


vsrl 


1  1  _  5 

Vmk  2  2  Vk  VA  ’ 


where  T  is  the  largest  value  of  y(x,y)  within  the  intersection  of  A  and  5. 


(3.6) 


The  standard  deviation  estimates  (3.5)  and  (3.6)  are  general  assuming  only  that  the 
range  of  v(x,y)  is  the  interval  [0,1]  and  that  the  sampling  points  are  random.  One  can 
obtain  other  estimates  by  also  taking  into  account  the  piecewise-constant  nature  of  the 
integrand  and  the  cellular  arrangement  of  the  sampling  points.  We  shall  derive  such 
estimates  in  later  sections  of  this  report.  However,  their  usefulness  is  likely  restricted  to 
simple  target  geometries,  because  the  application  of  the  formulas  require  an  analysis  of  the 
target  description. 
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4.  Cell  Size  Commensurate  with  Target  Description 

In  Section  2  we  introduced  a  subdivision  of  the  target  silhouette  in  nodal  areas  F}  with 
approximately  constant  v(x,y)  in  each  area.  The  proper  cell  size  for  the  numerical 
evaluation  of  vulnerability  integrals  depends  on  the  size  and  shape  of  these  areas.  On  one 
hand,  the  cells  should  be  sufficiently  small  to  capture  complicated  boundaries  of  the  nodal 
areas.  On  the  other  hand,  small  slender  extensions  from  a  nodal  area  do  not  contribute 
significantly  to  the  values  of  the  integrals.  We  characterize  the  smoothness  of  the 
boundaries  of  nodal  areas  by  a  parameter  that  expresses  the  fineness  of  the  details  of  the 
boundaries.  We  loosely  define  the  parameter  as  the  average  diameter  of  the  smallest 
relevant  convex  part  of  the  nodal  area  Ft  and  call  it  the  fineness  fj  of  the  nodal  area  Fr 
(The  relevance  of  parts  of  nodal  areas  are  discussed  at  the  end  of  Section  2).  The  fineness 
of  the  target  description  is  defined  as  the  smallest  f} : 

/=TU}  .  (4.1) 


In  practical  applications,  the  fineness  measures  f}  and  /  can  be  estimated  by  an  inspection  of 
the  projection  of  the  target  on  the  reference  plane  E.  The  fineness  measure  depends  on  the 
aspect  of  the  target,  on  the  type  of  the  threat  and  on  the  type  of  the  vulnerability  model.  It 
is,  therefore,  not  a  constant  for  a  given  three-dimensional  target  geometry  description. 
However,  the  order  of  magnitude  of  /  should  be  a  valid  indicator  for  the  smallest  diameters 
of  relevant  details  in  a  target.  An  example  of  /  is  discussed  in  Section  7. 

The  concept  of  the  fineness  measure  is  important  for  the  representation  of  target 
components  in  the  vulnerability  integral  approximations  (3.3)  and  (3.4).  If  the  distances 
between  the  sampling  points  (xu,yu)  are  larger  than  the  fineness  measure  then  relevant 
components  of  the  target  description  might  not  be  represented  in  the  integrals.  Also,  if  a 
cellular  sampling  point  arrangement  systematically  excludes  a  set  of  small  components  then 
the  standard  deviation  estimates  (3.5)  and  (3.6)  are  not  valid.  Therefore,  it  is  reasonable  to 
require  that  the  average  distance  between  sampling  points  is  smaller  than  a  fraction  (for 
instance,  smaller  then  one  half)  of  the  fineness  of  the  target  description.  On  the  other 
hand,  there  is  no  advantage  to  have  a  computing  grid  with  average  distances  between 
sampling  points  much  smaller  than  the  fineness  measure. 


An  average  distance  between  uniformly  distributed  sampling  points  in  a  plane  can  be 
defined  as  follows.  Let  k  sampling  points  be  located  in  an  area  At.  We  assign  to  each  point 
an  area  Ajk,  represent  this  area  by  a  circle  with  the  diameter  6,  and  define  8  as  the 
average  distance  between  the  points.  Explicitly,  for  given  sampling  area  Ak  and  number  k 
of  sampling  points  the  average  distance  is 


8  = 


(4-2) 


For  instance,  if  the  computing  cells  are  squares  with  the  side  length  s,  and  the  number  of 
samples  per  cell  is  k  then  the  average  distance  between  sampling  points  is 
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(4.3) 


8  =  ^  V4/ (it  k)  =  1.13  j/Vfc 

If  the  grid  is  triangular  with  the  side  length  then  the  corresponding  formula  is 

8  =  .v4  V(2  sin  60°)/(ir  k )  =  0.74  (Vk  .  (4.4)  I 

A  reasonable  value  of  the  average  distance  between  sampling  points  is 

S  *  jf  ■  (4.5) 


This  value  might  be  too  large  if  all  nodal  areas  are  square-like  or  circle-like.  On  the  other 
hand,  real  target  descriptions  characteristically  have  complicated  component  projections  and 
the  indicated  8  should  be  about  right  for  the  capturing  of  relevant  slender  nodal  areas. 


For  a  square  grid,  eq.  (4.5)  can  be  replaced  by 

s  <  0.44/  Vk 


or  by 


k  >  5.1 


(4.6) 

(4.7) 


Eq.  (4.6)  and  (4.7)  are  conditions  that  the  size  of  square  cells  and  number  of  sampling 
points  in  each  cell  are  commensurate  with  the  fineness  of  the  target  description. 

If  the  integrals  are  computed  by  a  general  Monte  Carlo  algorithm  without  a  grid 
pattern  then  the  condition  (4.5)  in  terms  of  the  number  k  of  samples  is 

* s  v  7  - 51 j  ■  <4-8> 

where  5  [m2]  is  the  sampling  area,  i.e.  the  silhouette  of  the  target  or  the  averaging  area  A. 


5.  Integration  Accuracy  Based  on  Nodal  Areas. 

In  this  section  we  estimate  the  standard  deviation  of  the  discrete  approximations  of 
vulnerability  integrals  for  the  case  where  the  integrand  can  be  represented  in  terms  of  nodal 
values  and  nodal  areas.  First  we  consider  the  vulnerable  area  Vs.  By  definition, 


;'“1  ' 


(5.1) 


where  at  is  the  fraction  of  the  nk  sampling  points  that  are  inside  the  nodal  area  Fr  If  the 
number  of  sampling  points  is  sufficiently  large  then  estimates  the  probability  that  a 
sampling  point  falls  within  the  nodal  area  FJt  that  is 


a 


(5.2) 
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and  the  standard  deviation  of  a,  is  (Ref.  1,  p.  168) 

\  1 1/2 

i  i  t* .  i  r .  I  I 

a(aj)  = 

Hence  the  standard  deviation  of  Vs  is 


F, 

Mil 

nk 

5 

s  J 

0(v’>  “  vh" 


s 

im  i 


V, 


5  J  S 


1  1 
Vnk  2 


M 


(5.3) 


1/2 


(5.4) 


Estimates  of  bounds  for  cr(Vj)  can  be  obtained  by  manipulations  outlined  in  the  Appendix. 
The  results  are 


/t 7  \  <•  5"  Vc 

<,(v’)  s  WT 


'V./itux  ^y  I 


V,  _ 

^y  ^FJ  min 


1/2 


5  Vc (  _  %  <  5  Vc 

2  V^k  2  ^  Vfjm“  Vf;rain  J  “  2  Vrt.  2 


(5.5) 


Even  the  least  accurate  bound  in  eq.  (5.5)  is  for  c< 4  smaller  than  the  general  bound  (3.5). 
If  rss 4  then  eq.  (5.4)  should  produce  a  better  bound  than  eq.  (3.5). 

For  the  average  vulnerability  we  have 

V.  -  V,  -  (M„  -  V,  )  =  -  V,  )  ,  (5.6) 

;-0  y-0  * 

where  b,  is  the  fraction  of  sampling  points  that  are  in  (FA),  =  F,C\  A,  and  ( FA)j/A  estimates 
the  probability  that  a  random  point  falls  within  ( FA), .  The  standard  deviation  of  the 
estimate  (FA),  /A  is 


(FA), 


1 

Vmk 


(FA),  (  _  (/vt),  ' 

A  {  A  , 


1/2 


or,  because  n/m  =  S/A 


(FA),  _  l  Vs 

,  A  Vnk  Va 


(FA), 

A 

A  J 

1/2 


(5.7) 


(5.8) 


Hence 


;=0 

Bounds  for  the  standard  deviation  are  in  analogy  to  eq.  (5.5) 


1/2 


(5.9) 
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~  ~Vnk  ~VA  ~i~~  [(^V^;i'«_V-'^VA-('>,V)Fjmin)]1  2  S 

-  2  VnF  W  [("V V) ("y V)f/min]  S 

1  s  Ve+1 
2Vfc  VX  2 


(5.10) 


The  accuracy  estimates  in  this  section  are  based  on  a  piecewise-constant  nature  of  the 
integrand  and  assume  a  random  distribution  of  sampling  points.  We  shall  see  in  the  next 
section  that  one  can  obtain  smaller  estimates  if  also  the  cellular  sampling  point  arrangement 
is  taken  into  account.  Therefore,  if  such  an  arrangement  is  used  and  the  cell  size  is 
commensurate  with  the  target  detail  then  eqs.  (5.4)  and  (5.9)  are  upper  bounds  of  standard 
deviations. 


6.  Integration  Accuracy  Based  on  Cells 

In  this  section,  we  estimate  the  accuracy  of  the  discrete  approximations  of  the 
vulnerability  integrals  in  the  case  where  the  sampling  points  are  arranged  in  a  cell  pattern 
and  the  cells  are  commensurate  with  the  target  description. 

Let  v,  be  the  average  value  of  v(x,y)  in  the  computing  cell  i.  Then  the  exact  value  of 
the  vulnerable  area  integral  (2.1)  is 

Vs  =  C  .  (6.1) 

i  =  l 


The  difference  between  the  approximation  (3.3)  and  Vs  is 


V 


s 


(6.2) 


The  summation  is  over  all  those  cells  which  intersect  with  the  silhouette  5.  Next,  we 
subdivide  the  n  computing  cells  into  two  groups.  Let  one  group  consist  of  7  cells  that  do  not 
contain  nodal  area  boundaries,  and  let  the  set  of  the  corresponding  cell  indexes  be  {/}.  Let 
the  parameters  of  the  other  group  (cells  containing  boundaries)  be  denoted  by  B  and  {b}, 
respectively.  Then  eq.  (6.2)  can  be  reformulated  as  follows: 


Vs-V's  =  C  2 

1  k 

1  k 

7  2va  - 

<<{<} 

*  /  =  ! 

k  /- 1 

The  first  sum  in  eq  (6.3)  is  close  to  zero  because  by  definition  the  effect  function  v(x,v)  is 
about  constant  within  each  cell.  We  assume  in  the  following  that  the  first  sum  can  be 
neglected. 
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The  second  sum  in  eq.  (6.3)  is 


S„=C 


(6.4) 


The  summation  is  over  the  5  cells  which  contain  component  boundaries.  Let  Av,  be  the 
range  of  v(x,y)  in  the  cell  i.  Then  j St  \  can  be  bounded  by 

|S*  |  SC  2  Av,  =  ^  £Av,  •  (6-5) 

i<{i}  «({b} 

A  cellular  arrangement  of  the  sampling  points  has  advantages  only  if  the  cell  size  s  is  small 
compared  to  the  target  description  fineness  measure  /.  We  therefore  assume  that 
commensurability  is  satisfied  by  a  large  margin: 

s  «f  .  (6.6) 

We  now  split  the  sum  (6.5)  over  all  boundary  cells  in  partial  sums  over  cells  containing  the 
boundary  of  each  nodal  area  Fj,  and  over  the  outer  boundary  of  the  integration  area  S,  and 
estimate  the  bound  of  |Sj  by  separately  bounding  each  of  these  partial  sums.  Let  pFj  be  the 
length  of  the  perimeter  of  Fj,  ps  be  the  length  of  the  outer  boundary,  Av^  be  the 
maximum  of  the  absolute  value  of  the  discontinuity  of  v(x,y)  across  the  boundary  of  F., 
hvpS  be  the  maximum  absolute  value  of  the  discontinuity  across  the  outer  boundary,  and 
Av„  be  the  maximum  absolute  value  of  all  discontinuities: 

Av  =  max  {  Av, }  .  (6.7) 

J 


Let  be  the  total  length  of  all  boundaries  within  the  integration  area: 


Ps 


U- 1 


(6.8) 


The  discretization  error  (6.3)  of  the  vulnerable  area  is  equal  to  the  absolute  value  of  5*  and 
can  be  bounded  by 


|s  Av,  <  f 

u{h} 


2  Pn  AV, 
U= 1 


+  Ps  *VpS 


(6.9) 


and  furthermore  by 

I  Vs~vs  I  s  *  Pz  Av,  s  px 

An  estimate  for  the  perimeter  pFj  is 


Pf,~  2  J- 

J  j 


1  +  i 


(6.10) 


(6.11) 
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The  formula  (6.11)  is  about  right  for  slender  F}  and  might  be  used  if  the  perimeter  pFJ  is 
not  available.  However,  in  general  estimates  of  the  pF)  can  be  easily  determined  by  an 
inspection  of  the  target  description. 


In  a  single  boundary  cell,  an  estimated  bound  for  the  standard  deviation  of  the 
discretization  is  (see  Appendix) 


a(S“)  *  C  Vk  2  =  Vk  2  4v‘  ' 

The  corresponding  estimate  of  the  standard  deviations  of  Sb  and  Vs  is 


-(Vj)  =  o(S„)  ^  2  (Av/)2 

(if  {b} 


M/2 


or 


,3/2 


-(v5)  = 


S Pfj  (AvPfj)2  +  Ps  (AvPs)2 

U-1 


1/2 


2  Vk 


Av_ 


(6.12) 


(6.13) 


(6.14) 


Eqs.  (6.13)  and  (6.14)  are  alternatives  to  eqs.  (5.4)  and  (5.5).  In  order  to  compare  both 
sets  of  formulas  we  use  the  relation  s~VS/n  and  obtain 


-(Vs) 


*3'2  (PxYn 
2  Vk 


Avp 


-i  .  gWjAv 

„y*  Vrt  2  ' 


(6.15) 


Thus,  for  sufficiently  large  n  the  estimates  (6.13)  or  (6.14)  are  smaller  than  the 
estimates  (5.4)  or  (5.5)  that  do  not  take  into  account  the  cellular  arrangement  of  sampling 
points. 

We  now  consider  the  average  vulnerability  integral  (2.5).  Let  its  exact  value  be 


V 


A 


(6.16) 


where  the  summation  is  over  the  m  cells  that  intersect  with  A.  The  difference  between  the 
discretization  (3.4)  and  V A  is 


^-v,-£i 

f  > 

T 2)  ('**')*  "  (”YV)/ 

1  " 
-is 

l| 

1 

& 

-W 

i-t 

/«1 

m  i-i 

1  *  “i  J 

Within  each  computing  cell  the  difference  is  bounded  by 

I  VA  ~  va  I  *  s  +  7,  Av,  )  ,  (6.18) 

m 

where  v,  and  y,  are  average  values  of  v(jc,y)  and  y(jt,y),  respectively,  in  the  cell  i,  and  Av( 
and  Ay,  are  the  corresponding  variations.  Let  {A}  be  the  set  of  cell  indexes  of  cells  that 
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(6.19) 


intersect  with  A.  The  first  term  in  eq.  (6.18)  provides  the  sum 

r  m  _  i  m  _ 

=  -j  2  ^  =  —  2  v<  A^ s  k A^i  • 

A  iml  m  jml 


The  second  term  in  eq.  (6.18)  contributes  nothing  if  the  cell  is  inside  a  nodal  area.  For  cells 
containing  nodal  boundaries,  we  can  estimate  the  total  contribution  to  the  difference  VA—VA 
by  using  the  same  steps  as  for  the  accuracy  estimates  of  Vs.  The  only  important  difference 
is  that  all  perimeter  lengths  and  maxima  pertain  to  the  averaging  area  A  instead  of  S.  We 
indicate  this  by  placing  a  tilde  over  the  corresponding  symbols.  The  result  is 


C  m  -  ,  * 

dia  ~  T  Av<  s  2 ~a  2  Pn  'tp'i  AV/  +  ip*  AV  ' 

1-1  [jrn 0 

where  ypFJ  and  ypA  are  the  largest  values  of  y^.y)  on  the  boundaries  pFJ 
respectively.  Eq.  (6.20)  can  be  bounded  further  by 


(6.20) 

and  pA, 


D 


2  A 


T  P*  yp  Ai>  s  7  W* 


(6.21) 


where  p^  is  the  total  length  of  boundaries  within  A  and  ip  is  the  largest  value  of  ■y(x,y)  on 
all  boundaries.  Combining  the  estimates  0U  and  D2A  one  obtains,  for  instance  the  bound 


i  V.-Vjs  max  [v,  A^]  + 
ie{A} 


(6.22) 


or  some  similar  expression  based  on  some  of  the  forms  of  0M  and  DlA .  In  order  to  expose 
the  trend  of  the  estimate,  we  introduce  the  quantity  y’  by  the  definition 


(6.23) 


For  linear  -y(x,y)  it  approximates  the  absolute  value  of  the  gradient  of  y  in  the  cell  i.  Then 
eq.  (6.22)  can  be  reformulated  as  follows 


max  [v,  *»;] 

|<{A} 


PXA 

A 


Ip  Avp  )  - 


(6.24) 


which  shows  that  the  estimate  of  the  difference  is  proportional  to  the  cell  size  s. 


To  obtain  the  standard  deviation  of  V A  we  follow  the  same  steps  as  for  the  derivation 
of  eqs.  (6.12)  through  (6.15).  The  result  can  be  formulated  by 


o(yA) 


1/2 


(6.25) 


or  in  one  of  the  following  equivalent  forms 
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^  =  vTvT 


max  [v,  A7J2  + 
a{A] 


s  Pu 
4  A 


(yP  *Vp)2 


°(VA)  * 


1 

rtv 4  Vnk 


(6.26) 

(6.27) 


The  estimate  (6.27)  is  smaller  than  (5.10)  if  n  is  sufficiently  large. 

The  estimates  in  this  section  presuppose  small  cell  sizes  (s«/)  whereas  the  estimates 
in  Section  5  depend  only  on  the  number  of  samples.  If  the  sample  number  n  is  increased 
then  eventually  the  estimates  in  this  section  become  better  (smaller)  than  those  of  Section  5. 
If  n  is  not  very  large  then  the  estimates  of  Section  5  might  be  smaller.  If  there  is  a 
difference  then  in  general  one  should  use  the  smaller  estimate.  The  estimates  at  the  end  of 
Section  3  are  most  conservative,  because  they  do  not  postulate  any  properties  of  the 
integrand  or  sampling  point  arrangement.  On  the  other  hand,  they  are  very  simple  to 
calculate  and  might  be  the  only  practical  formulas  if  the  target  description  is  complicated. 


7.  Application  Example 

In  the  previous  sections  we  have  derived  a  number  of  accuracy  estimates  and  criterions 
for  a  rational  choice  of  the  numerical  integration  method  and  of  the  level  of  detail  in  target 
description.  In  this  section  we  illustrate  the  application  of  these  results  with  a  simple 
artificial  example.  We  have  chosen  an  artificial  example  to  have  a  clear  presentation  of  the 
principal  aspects  of  the  criterions.  In  a  real  life  example  one  would  have  to  discuss  the 
possibly  complicated  details  of  the  target.  Such  discussions  are  necessarily  lengthy  and  do 
not  contribute  to  the  understanding  of  the  method. 

Table  1.  Sizes  and  Vulnerabilities  of  Nodal  Areas. 


n 

F  [m2] 

P  [m] 

/H 

V 
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m 
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0.05 

n 

.  ■  ■  1 1 1 
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0.8 

IE9 

0.0 

1.10 
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1.0 

B9 

0.05 
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5.33 
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0.15 
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0.05 

5 

1.00 

40.0 

0.05 

1.0 

0.0 

na 

2.57 

1 

BIB 

1.0 

0.0 

El 

4.33 

my 

m 

0.4 

0.05 

5  =  A  =  11  m2,  ps  —  15  m, 


p±  —  27.34  m,  f  ~  fi  -  0.15  m, 
plt  =  67.34  m,  /t  =  f2i  =  0.05  m. 

Figure  2  shows  the  projection  of  a  target  description  onto  the  reference  plane.  The 
corresponding  numerical  data  are  listed  in  Table  1.  We  shall  use  these  data  to  compute  the 


-  17  - 


average  vulnerability  whereby  the  averaging  is  over  the  whole  silhouette.  (Hence,  A-S 
and  -y«*l).  The  first  four  entries  in  the  table  contain  the  parameters  of  the  four  nodal  areas 
shown  in  Figure  2.  The  last  three  entries  we  shall  discuss  later.  Using  the  data  for 
7=1,. ..,4  we  compute  with  eqs.  (2.13)  and  (2.14)  the  average  vulnerability  V  A  of  the  target 
consisting  of  the  first  four  components 

VA  =  0.430  ,  (7.1) 

and  the  corresponding  intrinsic  spread  DA  of  V A: 

Da  =  0.043  .  (7.2) 


The  commensurate  average  distance  between  sampling  points  equals  0.075  m 
according  to  eq.  (4.5).  If  a  square  grid  is  used  with  one  sampling  point  per  cell  then  the 
side  length  of  the  cell  should  be  less  than  0.06  m  (see  eq.  (4.6)  ).  If  the  side  length  is 
0.1  m  then,  according  to  eq.  (4.7),  one  should  take  two  samples  from  each  cell.  However, 
the  criterions  based  on  the  target  fineness  measure  implicitly  assume  that  the  fineness  is 
typical  for  the  target  description  and  that  the  nodal  fineness  measures  are  diameters  of 
relevant  areas.  In  the  present  example,  the  fineness  is  determined  by  the  width  of  the 
isthmus  between  the  nodal  areas  F2  and  F3.  It  is  obvious  that  a  failure  to  represent  the 
isthmus  in  the  vulnerability  integral  would  not  significantly  affect  the  average  vulnerability. 
Therefore,  a  better  value  for  the  fineness  measure  is  /=0.5  m  which  approximately 
represents  typical  diameters  of  the  narrow  parts  of  the  nodal  areas  if  the  isthmus  between 
and  F3  is  neglected.  Then  the  commensurate  cell  size  according  to  eq.  (4.6)  should  be 
less  than  0.2  m,  that  is,  a  cell  size  of  0.1  m  is  adequate. 

We  test  the  relevance  of  the  isthmus  by  using  eq.  (2.15)  and  assuming  that  the  area  of 
the  isthmus  is  added  to  F2.  (Adding  to  F3  has  a  smaller  effect).  The  effect  of  this 
replacement  on  the  average  vulnerability  is  less  than 

T  -j  (  v2  ~  v4  )  =  (0.15  0.5)- -7-O.6  =  0.0041  ,  (7.3) 

A  i  1 

which  is  much  smaller  than  the  intrinsic  error  DA  of  V A. 


Let  the  cell  size  for  the  numerical  integration  be  r  =  0.1  m,  and  let  the  number  of 
samples  from  each  cell  be  k—1.  Then  the  estimate  (3.6)  for  the  standard  deviation  of  the 
numerically  computed  average  vulnerability  VA  from  the  true  value  VA  is 

"(Ya)  s  777  [(r  -  va)  Va]V2  =  -■—9: =  0.0172  (7.4) 

which  is  about  40%  of  the  intrinsic  error  D A .  More  accurate  are  the  estimates  (5.9)  and 
(6.26).  With  eq.  (5.9)  we  obtain 


1 

Vnk 


0.2845 

Vnk 


0.0858-.*  , 


(7  5) 
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where  we  have  used  the  relation  Vn  =  VS/s  =  vTl/j.  This  estimate  shows  that  an 
integration  with  the  cell  size  5  =  0.1  m  is  adequate  because  the  standard  deviation  of  the 
numerical  integration  is  only  about  20%  of  the  intrinsic  error  DA  of  the  average 
vulnerability,  eq.  (7.2).  Even  a  cell  size  5  =  0.2  m  is  sufficiently  accurate  in  this  example 
and  would  require  only  a  fourth  of  the  computing  time. 

The  first  term  in  the  standard  deviation  estimate  formula  (6.26)  equals  zero,  because 
in  our  example  y  is  constant  and  Ay(  =  0.  The  second  term  yields  for  5  =  0.1 

—  r^/2  ( 3/2  \/'J  7 

a(VA)  ~  WklA  '  Av'  =  -  Vi!  --0'75  =  0'1782'*3/2  =0-0056  .  (7.6) 


Thus,  the  more  accurate  formulas  yield  in  this  example  estimates  that  are  smaller  than  the 
general  estimate  (7.4)  by  a  factor  of  0.5  to  0.3. 

The  bound  of  the  error  of  the  numerical  integration  can  be  commuted  using  eq.  (6.22): 

i  VA-VA  |  s  -Ai,  =  ^f^-0.75  -  1.864-5  .  (7.7) 


The  more  accurate  eq.  (6.20)  produces 


V  -V'  i<  —:L~ 

A  A  1  2  A 


2 PrAvrn  +  Ps^vts 
1/-1 


=  ^--27.688  =  1.258-5 
22 


(7.8) 


for  5  =  0.1  we  obtain  a  bound  that  is  three  times  as  large  as  the  intrinsic  error  DA  of  V A. 
Therefore,  if  guaranteed  error  bounds  arc  important  then  the  cell  size  s  should  be  about 
0.03  m.  However,  because  the  best  estimate  of  the  standard  deviation  of  VA  is  only  about 
13%  of  the  intrinsic  error  DA,  the  large  estimate  of  the  bound  likely  is  too  pessimistic  and 
can  be  ignored. 

We  now  assume  that  the  target  contains  in  addition  to  the  first  four  components  a  very 
sensitive  fifth  component  with  a  small  diameter  (e  g.,  a  fuel  line),  and  determine  whether 
that  component  should  be  included  in  the  target  description.  Let  the  effective  diameter  of 
the  component  be  2  cm,  its  length  be  10  m,  and  let  it  be  located  in  the  component  area  4. 
The  presented  length  of  the  line  is  rather  large:  it  is  two  times  as  long  as  the  horizontal 
dimension  of  the  silhouette  and  about  three  times  as  long  as  the  length  of  the  area  j- 4. 
We  use  eq.  (2.15)  and  estimate  the  effect  of  the  additional  component  on  the  average 
vulnerability  The  presented  area  of  the  component  is  r5  =  0. 02  10=0.2  m2  and 

T,  ■  j  ■  (v,-v4)  =  0.2—  0. 6  =  0.011  .  (7.9) 

This  is  about  one  fourth  of  the  intrinsic  error.  Therefore,  the  component  is  not  relevant 
and  the  inclusion  of  such  a  component  in  the  target  description  is  a  waste  of  resources.  On 
the  other  hand,  if  such  a  component  is  already  present  in  the  target  description,  then  it  can 
be  ignored  in  the  sense  that  the  cell  size  need  not  be  made  commensurate  with  the  fineness 
measure  of  the  component. 
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To  have  an  example  where  a  change  of  the  target  description  is  indicated  we  assume 
that  the  additional  fuel  line  has  an  effective  diameter  of  5  cm  and  a  length  of  20  ra.  It  is 
listed  as  the  fifth  component  in  the  Table  1.  Its  effect  equals  0.055  and  is  about  the  same  as 
the  intrinsic  error,  eq.  (7.2).  (Also,  the  effective  presented  area  of  the  line  is  formidable:  it 
is  about  the  same  as  the  area  F3  or  about  10%  of  the  silhouette).  Because  the  fifth 
component  is  assumed  to  be  located  in  the  nodal  area  j=  4,  the  area  and  its  perimeter  must 
be  modified.  We  assume  for  simplicity  that  the  fineness  measure  /4  is  not  changed,  and  list 
the  other  new  parameters  of  F^  in  the  line  j=  4t  of  Table  1.  Because  vj=v2,  the  new 
component  belongs  to  the  nodal  area  F2  the  parameters  of  which  are  changed  and  listed  in 
the  line  j=  2t-  Finally,  the  total  length  of  nodal  area  boundaries  />s  is  increased  to  p£t  and 
the  new  fineness  measure  of  the  target  description  is  /t=/5=/2t.  The  new  table  of  the 
target  description  now  consists  of  the  entries  j=  1,  2f,  3  and  4t. 

Using  the  modified  table,  one  obtains  for  the  average  vulnerability 

VAt  =  0.485  .  (7.10) 


Hence,  by  adding  the  fuel  line  to  the  target  description  the  average  vulnerability  is 
increased  by  0.055.  The  intrinsic  spread  of  the  new  result  is 

DA,  =  0.038  .  (7.11) 


The  spread  D„t  is  smaller  than  DA,  because  the  added  component  has  an  error  free  effect 
function  value  v5**l. 


The  standard  deviation  of  VA1  is,  if  computed  with  eqs.  (3.6),  (5.9)  or  (6.26), 
respectively, 


v(vA1)  =  =  o.oi5l  . 

_/,7  ^  0.3093  0  0933  n  „nn, 

a(v-)  =  -  vsr  =  ~vT  '  -  0  009 J 


■t3/2-V67.34 

Vk  2  11  V 


afV^)  ss  ~  0.75  =  0.2798  =  0.0088  . 


(7.12) 

(7.13) 

(7.14) 


Thus,  the  cell  size  .v  =  0.1  m  is  adequate  for  the  modified  target  in  the  sense  that  the  ic  .u’it 
of  the  numerical  integration  likely  will  be  within  the  intrinsic  spread  of  the  average 
vulnerability,  and  that  the  standard  deviation  of  the  integral  is  less  than  209 i  of  the 
vulnerability  increase  0.055  due  to  the  new  component 


Next  we  compute  estimates  for  the  bound  of  the  difference  between  V A,  and 
Using  eq.  (6.22)  one  obtains 

I  VAf  -  VAi  I  s  -^y3-  0.75  -  4.591  r  .  (7.15) 


and,  with  the  more  accurate  eq.  (6.20), 
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(7.16) 


I  VAf  -VAt\s  3.440-i  . 

Therefore,  to  guarantee  a  numerical  integration  accuracy  better  than  the  effect  of  the  fifth 
component  one  needs  a  cell  size  of  about  0.055/3.440=0.016  m.  This  is  about  the  same  as 
the  commensurate  cell  size  which  is  0.020  m  according  to  eq.  (4.6).  For  the  cell  size 
0.02  m,  the  number  of  cells  is  about  25  times  larger  than  for  5  =  0.10  m.  We  note  in  passing 
that  sampling  of  25  random  points  per  cell  results  in  the  same  increase  of  computing  time 
and  a  substantial  reduction  of  the  standard  deviation  estimate  but  does  not  guarantee  that 
the  effect  of  the  fifth  component  is  correctly  computed,  because  the  random  sampling  points 
might  not  fall  within  the  projection  of  the  component. 

In  summary,  a  cell  size  of  0.1  m  with  one  random  sample  per  cell  is  adequate  for  the 
computation  of  the  average  vulnerability.  Likely  even  the  effect  of  the  (long  and  thick)  fuel 
line  (the  fifth  component)  will  be  correctly  determined.  However,  because  the 
commensurability  condition  is  not  satisfied  a  detection  of  the  effect  of  the  fifth  component 
is  not  guaranteed.  If  a  guaranteed  computation  of  that  effect  is  important  then  the  cell  size 
must  be  reduced  at  the  cost  of  a  corresponding  increase  of  the  computing  time  by  a  factor 
of  25  or  more. 

Finally,  the  example  illustrates  the  tediousness  of  the  application  of  the  theoretically 
better  formulas  for  error  estimates.  Usually,  the  simple  standard  deviation  estimate, 
eq.  (7.4)  or  (7.12),  will  suffice  to  establish  adequacy  of  the  numerical  integration,  and  a 
detailed  analysis  of  the  target  description  will  not  be  necessary.  The  significance  of  the 
more  accurate  formulas  is  mainly  theoretical,  showing  that  a  cellular  point  arrangement  is 
indeed  better  than  random  sampling. 

8.  Summary  and  Conclusions 

We  considered  in  this  report  the  computation  of  target  system  vulnerabilities  by 
numerical  integration  of  local  effect  functions  over  a  reference  plane.  Accuracy  is  affected 
by  three  sources  of  errors:  an  "intrinsic"  error  caused  by  inaccurate  models  of  component 
response,  an  error  caused  by  inadequate  target  geometry  description  and  a  numerical 
integration  error.  In  a  reasonable  algorithm  the  magnitudes  of  these  three  errors  should  be 
of  the  same  order,  because  the  final  accuracy  is  dominated  by  the  largest  error  and  a 
reduction  or  elimination  of  one  or  both  of  the  smaller  errors  does  not  change  the  order  of 
magnitude  of  the  total  error.  In  this  report,  we  assume  as  given  the  "intrinsic"  error  which 
is  determined  by  the  accuracies  of  component  vulnerabilities.  Errors  from  the  other  two 
sources  can  be  controlled  by  proper  choices  of  target  description  detail  and  numerical 
integration  procedures,  respectively.  Error  estimates  provided  in  this  report  help  one  to 
make  such  choices. 

We  assume  that  the  integration  is  done  with  a  Monte  Carlo  method  as  follows.  First, 
the  integration  area  is  overlaid  with  a  computing  net  consisting  of  square  cells  with  the  cell 
side  lengths  5  =  0.1  m.  Next,  the  integrand  is  computed  at  k  random  sampling  points  in  each 
computing  cell,  and  the  integral  is  obtained  by  weighted  adding  of  the  sampled  integrand 
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values.  The  usual  computation  is  done  with  Jfc  =  l  sampling  points  from  each  computing  cell. 

The  following  particular  aspects  might  be  treated  based  on  the  results  presented  in  this 
report: 

(a)  Estimates  and  bounds  of  the  overall  accuracy  of  vulnerability  calculations. 

(b)  Relevance  of  particular  components. 

(c)  Rational  choice  of  computing  mesh  based  either  on  the  fineness  of  the  target 
description  or  on  the  achievable  accuracy. 

Formulas  which  may  be  used  to  address  these  problems  are  listed  in  Table  2  by  application 
area. 

Table  2.  Numbers  of  Formulas  by  Application. 


Application 

Vulnerable 

area 

Average 

vulnerability 

Intrinsic  error 

(2.8),  (2.10) 

(2.12),  (2.14) 

Standard  deviation 
of  integral 

(3.5), 

(5.4),  (5.5), 

(6.13),  (6.14),  (6.15) 

(3.6), 

(5.9),  (5.10), 

(6.25),  (6.26),  (6.27) 

Bounds  of  integration 

error 

(6.5),  (6.9),  (6.10) 

(6.19),  (6.20), 

(6.22),  (6.24) 

Component  relevance 

(2.11) 

(2.15) 

Commensurate 
mesh  size 

(4.6) 

Three  general  conclusions  can  be  drawn  from  the  discussions  in  this  report.  First,  a 
statement  about  the  computing  cell  size.  The  standard  deviation  of  the  Monte  Carlo 
computation  of  the  vulnerability  averaged  over  the  silhouette  is  (see  eq.  (3.6)  ) 

Assuming  5  =  10  m2  as  typical  for  the  order  of  magnitude  of  an  armored  vehicle  silhouette, 
and  setting  j  =  0.1,  one  obtains 

c x(Vj  <  0.16-.V  =  0.016  . 

Hence,  for  a  cell  size  of  j  =  0.1  m  the  standard  deviation  of  the  integral  is  of  the  order  of 
the  second  digit  of  the  average  vulnerability.  This  is  likely  less  than  the  intrinsic  error  of 
the  average  vulnerability  VA  due  to  model  inaccuracies  and,  therefore,  a  cell  size  of  0.1  m  is 
adequate  for  the  computation  of  directional  average  vulnerability  of  a  typical  tank. 

The  second  conclusion  pertains  to  the  importance  of  detail  in  the  target  description. 
Minute  details  are  necessary  if  one  wants  to  predict  the  effect  of  a  single  shot  on  a 
particular  target.  For  any  average  vulnerability  which  is  defined  by  an  integral  over  an  area 
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in  the  reference  plane,  the  contribution  of  a  detail  is  directly  proportional  to  the  relative 
size  of  its  presented  area.  If  the  intrinsic  inaccuracy  of  the  average  vulnerability  indicator  is 
of  the  order  of  5%-10%,  a  typical  value,  then  any  detail  with  an  effective  presented  area 
of  less  than  10%  of  the  averaging  area  is  suspect  of  not  being  relevant.  (The  effective 
presented  area  is  defined  as  that  part  of  the  reference  plane  where  the  threat  has  an  effect 
on  the  component,  see  page  7).  Such  a  detail  might  be  removed  from  the  target  description 
without  penalty  as  shown  in  the  example  in  Section  7.  However,  the  removal  should  be 
judicious,  because  a  combination  of  several  irrelevant  components  might  well  be  relevant. 
On  the  other  hand,  if  one  wants  to  investigate  the  effect  of  a  small  detail  of  the  target  by 
comparing  the  values  of  the  corresponding  vulnerability  integrals,  then  the  mesh  size  must 
be  commensurate  with  the  detail.  (A  simple  increase  of  the  number  k  of  points  per  cell 
does  not  guarantee  that  the  effect  of  a  detail  is  correctly  computed).  This  can  increase  the 
computing  times  by  orders  of  magnitude  unless  the  algorithm  is  arranged  such  that  a  fine 
net  is  used  only  in  the  area  of  the  detail. 

Finally,  for  a  piecewise  constant  integrand,  cellular  arrangements  of  sampling  points 
theoretically  yield  more  accurate  results  than  random  sampling  if  the  ceil  size  is  much 
smaller  than  the  fineness  of  the  target  description. 
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LIST  OF  SYMBOLS 


A  -  area  of  averaging  in  the  reference  plane  E,  m2, 

{A}  -  set  of  cell  indexes  of  cells  intersecting  with  A, 

dj  -  fraction  of  sampling  points  in  the  nodal  area  FJt  dimensionless, 

{6}  -  set  of  cell  indexes  of  cells  containing  boundaries  of  nodal  areas, 

5  -  cardinality  of  {b},  dimensionless, 

C  -  cell  size  in  computing  grid,  m2, 

c  -  number  of  nodal  areas  within  the  silhouette  S,  dimensionless, 

DA  -  bound  of  bVA,  dimensionless, 

Ds  -  bound  of  8ys,  m 2, 

E  -  reference  plane  normal  to  the  direction  of  threat, 

e  -  number  of  nodal  areas  within  the  averaging  area  A,  dimensionless, 

Fj  -  nodal  area  j,  m2, 

(FA)j  -  intersection  of  Fj  with  A,  m2, 

/  -  fineness  of  the  target  description,  m, 

fj  -  fineness  of  the  description  of  the  nodal  area  F,,  m, 

£(•*>>0  -  density  of  averaging  weight,  m-2, 

k  -  number  of  samples  in  a  computing  cell,  dimensionless, 

m  -  number  of  cells  intersecting  with  the  averaging  area  A,  dimensionless, 

n  -  number  of  cells  intersecting  with  the  silhouette  5,  dimensionless, 

Px  -  total  length  of  boundaries  within  the  silhouette  S,  m, 

Pxa  ~  total  length  of  boundaries  within  the  averaging  area  A,  m, 

pA  -  perimeter  of  the  averaging  area  A,  m, 

Pfj  ~  perimeter  of  the  nodal  area  Fjt  ra, 

ps  -  perimeter  of  the  silhouette  S,  m, 

S  -  area  of  target  silhouette  in  the  reference  plane  £,  m2, 

s  -  side  length  of  a  square  cell  in  a  computing  grid,  m, 

Tj  -  effective  presented  area  of  a  target  component,  m2, 

{r}  -  set  of  cell  indexes  of  cells  completely  inside  a  nodal  area, 

1  -  cardinality  of  {/},  dimensionless, 

V A  -  average  vulnerability,  dimensionless, 

V A  -  approximate  average  vulnerability,  dimensionless, 

SVA  -  bound  of  error  of  VA  due  to  error  of  v(x,y),  dimensionless, 

Vj  -  vulnerable  area,  m2, 

V s  -  approximate  vulnerable  area,  m2, 

8VS  -  bound  of  error  of  Vs  due  to  error  of  v(x,y),  m2, 

v(*i>0  -  effect  of  a  threat  impinging  at  (x,y),  dimensionless, 

vit  -  sample  v(xu,yu),  dimensionless, 

vFJ  -  nodal  (average)  value  of  v(x,y)  in  the  nodal  area  Fjf  dimensionless, 

v T,  -  average  value  of  v(x,y)  in  the  effective  area  Tp  dimensionless, 
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v,  -  average  value  of  v(x,y)  in  cell  i,  dimensionless, 

v0J  -  average  value  of  v(x,y)  in  T}  without  the  component,  dimensionless, 

Av,  -  discontinuity  of  v(x,y)  in  cell  /,  dimensionless, 

Avf  -  maximum  discontinuity  of  v(x,y),  dimensionless, 

Av,  -  maximum  discontinuity  of  v(x,y)  in  A,  dimensionless, 

Av^  -  maximum  discontinuity  of  v(x,y)  on  pA,  dimensionless, 

A vpFJ  -  maximum  discontinuity  of  v(x,y)  on  pF],  dimensionless, 

&vpFJ  -  maximum  discontinuity  of  v(x,y )  on  pFJ  in  A,  dimensionless, 

Av^  -  maximum  discontinuity  of  v(x,y)  on  ps,  dimensionless, 

8v(x,y)  -  error  bound  of  v(jr,y),  dimensionless, 

8vf/nux  -  maximum  value  of  8v(x,y)  in  Ft,  dimensionless, 

-  maximum  value  of  8v(x,y)  in  7y,  dimensionless, 
x,y  -  Cartesian  coordinates  in  the  reference  plane  E ,  m, 

I"  -  largest  value  of  y(x,y)  within  AQ  S,  dimensionless, 

y(x,y)  -  normalized  density  of  averaging  weight,  dimensionless, 
yp  -  maximum  of  y(x,y)  on  nodal  boundaries,  dimensionless, 
y0  -  sample  y(x lj,ylj),  dimensionless, 

yp  -  maximum  y(x,y)  on  the  outer  boundary,  dimensionless, 

ypA  -  maximum  y(x,y)  on  the  perimeter  pA,  dimensionless, 

ypF]  -  maximum  y(x,y)  on  the  perimeter  pFjt  dimensionless, 

ypS  -  maximum  y(x,y)  on  the  perimeter  ps,  dimensionless, 

A-y,  -  absolute  variation  of  y(jt,y)  in  cell  i,  dimensionless, 

(yv)l  -  average  value  of  y(x,y)v(x,y)  in  cell  i,  dimensionless, 

(yv)FJ  ~  average  value  of  y(x,y)v(x,y)  in  the  nodal  area  Fj,  dimensionless, 

( yv)TJ  -  average  value  of  'Y(x,y)v(x,y)  in  the  component  area  Tjt  dimensionless, 

(7v)oy  ~  average  value  of  (-y  v)  in  area  T p  without  the  component,  dimensionless, 
-  maximum  value  of  [-y8v]  in  Fj,  dimensionless, 
b8  v]r;mM  -  maximum  value  of  [-y8v]  in  Tjt  dimensionless, 
yj  -approximate  |  grad -y  |  in  the  cell  i,  1/m, 

8  -  average  distance  between  sampling  points,  m, 

o-^)  -  estimated  standard  deviation  of  VA ,  dimensionless, 

a(Vs)  -  estimated  standard  deviation  of  Vs,  m2. 
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Appendix.  Estimate  of  Monte  Carlo  Standard  Deviation. 

Let  h(x)  be  integrable  over  the  x-interval  [a,b],  and  let  it  have  the  range 


*  “  T^7  /*<*>  dx  ■ 


We  estimate  h  by  chosing  k  random  points  Xjt[a,b]  and  computing  the  sample  mean 

=  t  £  • 


(A.l) 


(A.  2) 


The  expected  value  of  Hk  is  h  and  an  estimate  of  the  variance  of  the  sample  mean  Hk  is 
(Ref.  1,  p.  147) 


Var (Hk)  =  ~  2  Var(h(xj))  =  ±  J,(h (xj)  -  Hk )2  = 
k  j=  1  *  ;=1 


=  F  2*’W  -  7»!  - 

y-i 
*  y=l 

A  simple  upper  bound  of  Var(//t)  is  obtained  by  taking  the  largest  term  in  the  sum: 

Var (Hk)  <  j  max  [h(Xj)-Hk]2  <  j  max  {  {hm„-Hk)\  (Hk-hmin?  }  . 

A  generally  more  accurate  bound  is 


(A.  3) 


(A. 4) 


-  {Hk-KJ  -  jWt-KJ2  =  (A.5) 

k  j- 1 

To  show  that  V  is  an  upper  bound  of  Var (Hk)  we  compute  the  difference 
D  =  Var (Ht)  -  U  =  ±  £  [(A^)- hmin)2  -  (h(x})-hmin)(hmtt-hmin)]  .  (A.6) 

k  j-\ 

Eq.  (A.6)  shows  that  D  is  in  genera]  negative  and  vanishes  only  if  the  range  of  h(x)  is 
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restricted  to  the  set  {hmiB,hmtt}.  Therefore, 

Var (Hk)  s  1  (hm„-Hk)(Hk-hmia)  =s  1  i  (hmt-hmiay  ,  (A.7) 

since  Ht  can  only  have  values  within  the  interval 

The  standard  deviation  of  a  Monte  Carlo  calculation  thus  can  be  estimated  using 
eq.  (A. 4)  by 

a(Hk)  S  ^  raaX  {  Kxx~Hk<  Hk~hml«  1  S  -^(hmxx~hm\a)  »  (A-8) 

or,  using  the  more  accurate  bound  (A.7)  by 

CT(fy)  S  ( (^mtx (Hk  “  ^min) ) 1  2  ~  y  (hmxx~hm\n)  •  (A9) 

These  formulas  are  not  restricted  to  integration  in  one  dimension,  because  we  have  only 
assumed  integrability  and  used  the  maximum  and  minimum  properties  of  h(x).  Therefore, 
the  same  formulas  apply  to  the  calculation  of  integrals  in  any  dimension,  in  particular  to  the 
surface  integrals  in  vulnerability  analysis. 
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x  — Axis 


Figure  1.  Integration,  in  1—  D 
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Figure  2.  Target  projection  on  Fhe  reference  plane 
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2  Commandant 

US  Army  Field  Artillery  Center  &  School 
ATTN:  ATSF-CO-MW,  B.  Willis 
Fort  Sill,  OK  73503-5600 


1  Office  of  Naval  Research 

ATTN:  Code  473,  R.  S.  Miller 
800  N.  Quincy  Street 
Arlington,  VA  22217-9999 


3  Commandant 

US  Army  Armor  School 

ATTN:  ATZK-CD-MS,  M.  Falkovitch 

Armor  Agency 

Fort  Knox,  KY  40121-5215 

2  Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA  62R 
SEA  64 

Washington,  DC  20362-5101 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-954-Technical  Library 
Washington,  DC  20360 

1  Naval  Research  Laboratory 

Technical  Library 
Washington,  DC  20375 

1  Commander 

Naval  Surface  Warfare  Center 
ATTN:  Code  DX-21  Technical  Library 
Dahlgren,  VA  22448-5000 
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1  Commander 

Naval  Weapons  Center 

ATTN:  Information  Science  Division 
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2  Superintendent 

Naval  Postgraduate  School 
Department  of  Mechanical  Engineering 
Monterey,  CA  93943-5100 

1  Commander 
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1  AFSC/SDOA 

Andrews  AFB,  MD  20334 


33 


No.  of 

Copies  Organization 

1  AF  Astronautics  Laboratory 
AFAL/TSTL,  Technical  Library 
Edwards  AFB,  CA  93523-5000 

1  AFATL/DLYV 

Eglin  AFB,  FL  32542-5000 

1  AFATL/DLXP 

Eglin  AFB,  FL  32542-5000 

1  AFATL/DUE 

Eglin  AFB,  FL  32542-5000 

1  NASA/Lyndon  B.  Johnson  Space  Center 

ATTN:  NHS22  Library  Section 
Houston,  TX  77054 

1  FTD/NIIS 

Wrighl-Patlcrson  AFB,  OH  45433 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Technical  Information 
Department  L-3 
P.  O.  Box  808 
Livermore,  CA  94550 

2  Director 

Los  Alamos  Scientific  Laboratory 
ATTN:  Document  Control  for  Reports 
Library 

P.  O.  Box  1663 

Los  Alamos,  NM  87544 

1  Director 

Sandia  Laboratories 
ATTN:  Document  Control  for  3141 
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1  Director 
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Administration 

Scientific  &  Technical  Information 
Facility 

P.  O.  Box  8757 
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1  Aerospace  Corporation 
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1  The  Boeing  Company 

ATTN:  Aerospace  Library 
P.  O.  Box  3707 
Seattle,  WA  98124 

1  Director 

Institute  for  Defense  Analyses 
ATTN:  Library 
1801  Beauregard  Street 
Alexandria,  VA  22311 

1  Battelle  Memorial  Institute 

ATTN:  Technical  Library 
505  King  Avenue 
Columbus,  OH  43201-2693 

1  Johns  Hopkins  University 

Applied  Physics  Laboratory 
ATTN:  Jonathan  Russ 
John  Hopkins  Road 
Laurel,  MD  20707-0690 

1  Pennsylvania  State  University 

Department  of  Mechanical  Engineering 
ATTN:  K.  Kuo 

University  Park,  PA  16802-7501 

1  SRI  International 

Propulsion  Sciences  Division 
ATTN:  Technical  Library 
333  Ravcnswood  Avenue 
Menlo  Park,  CA  94025-3493 

I  Rensselaer  Polytechnic  Institute 

Department  of  Mathematics 
Troy,  NJ  12181 

1  Eichelbcrgcr  Consulting  Company 

ATTN:  Dr.  R.  Eichelbergcr,  President 
409  West  Catherine  Street 
Bel  Air,  MD  21014 

1  FMC  Corporation 

Aberdeen  Regional  Office 
ATOM:  Lcland  Watcrmcicr 
Bel  Air,  MD  21014 
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M.  Smith 
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ATTN:  AMSTE-SI-F 
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1  Cdr,  CSTA 
ATTN:  STECS 


35 


Imp.ntionai.ly  u. 


n  plank. 


36 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number _ brl-tr-3084 _ _  Date  of  Report  APR  _ 

2.  Date  Report  Received _ 


3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 
for  which  the  report  will  be  used.) _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  Genterifl  Comments.  What  do  you  thrnk  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Name 

CURRENT  Organization 

ADDRESS _ 

Address 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 

OLD 

ADDRESS 

Organization 

Address 

City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


FOLD  HERE 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN.  SLCBR-DDT 

Abertlcen  Proving  Ground,  MI)  2  IDs  50A6 

OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CUSS  PERMIT  No  0001,  APG.  MD 


POSTAGE  WILL  B£  PAID  BY  ADDRESSEE 


Director 

U.S.  Amiy  Ballistic  Research  Laboratory 

ATl’N:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-9989 


